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Abstract 

Some thermodynamical magnitudes of interest in a pure neutron plasma are studied within 
the framework of the non-relativistic Brueckner-Hartree-Fock approximation at finite density and 
temperature. We use Skyrme and Gogny forces to describe such a neutron plasma and study the 
main differences that arise in these two effective parametrizations of the nuclear interaction when 
a strong magnetic field induces a permanent magnetization in the gas. The existence of a non-zero 
permanent spin polarization in a neutron plasma is explored in the density-temperature parameter 
space. We find that for moderate temperatures and in the low density range up to densities ~ 0.5po 
both parametrizations predict that as density decreases an increasingly strong magnetization is 
allowed. In the range 0.5/9o p ^ 3/3o there is an approximately constant polarization that can 
be as big as ~ 12% for the maximum allowed interior magnetic field B ~ 10 18 G. For higher 
densities there is a dramatic difference in the polarization trend followed by Skyrme an Gogny 
forces. While the former predict a ferromagnetic phase transition, the Gogny forces prevent it 
keeping the magnetization below 5%. 
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I. INTRODUCTION 



The presence of magnetic fields in neutron stars and their effect in the thermodynamics 
of plasmas in the inner regions of neutron stars and in proto-neutron stars has been a 
subject of study for long time, see jl| and references therein. Since the pioneering work of 
Woltjer [2] who predicted very intense magnetic fields in neutron stars as a result of flux 
conservation from the projenitor star until the recent experimental measurements of pulsar 
spin-down rates or X-ray spectral lines from protons inmersed in such intense magnetic 
fields [3| much work has been done. Experimental insight into the field has been obtained 
from meaurements ranging from radio pulsars where magnetic fields are believed to be in 
the range B as 10 9 — 10 12 G to the so called "magnetars" predicted to exist by Duncan 
and Thompson in 1992 |4| and first detected in 1998 [5] where B magnetar w 10 14 - 10 15 G. 
These field strengths strongly decay in a relatively rapid time of 10 3 — 10 4 years. Two types 
of astronomical objects seem to be golden candidates for magnetars: soft gamma repeaters 
(SGRs) and anomalous X-ray pulsars (AXPs). Experimental observations confirm the idea 
that both are highly magnetized young neutron stars and can be linked to their original 
Supernova remnants [6]. Magnetars are believed to be around 10% of the whole population 
of neutron stars so they are relatively frequent objects although only a handful of both 
have been observed [l|. Other unconfirmed indications of strong magnetic fields include 
the possibility of an initial kick-off velocity due to an anisotropic neutrino-pair emission as 
suggested in |7J, |8J] . 

From a theoretical point of view, to study how matter behaves in the presence of such 
strong magnetic fields some works have dealt with the possibility that nuclear plasmas may 
be partially magnetized 
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It is widely known that Skyrme models 



have a characteristic pathology allowing a ferromagnetic transition at high densities even 



151 ] . However, this transition is not confirmed when 



in the absence of magnetic fields 

—ions are penned nsin g realise interactions Q. In the „nte r sne.ls of nentron 
stars one should focus in the low density limit of nuclear plasmas where clustered phases 



may develop 



17 



18 



19 



201 ] . For a magnetized system thus, we may have microscopic 



spin saturated spatial structures similar to magnetic domains. If present, they may largely 
enhance the axial response in the neutrino opacities relevant to proto-neutron star early 
cooling. Since neutrino transport seems to be a crucial ingredient in the simulation of 
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supernova core collapse it is important to size the effect of the existence of a strong magnetic 
field in the microphysics of these plasmas. 

The present work performs a comparative study of thermodynamical observables using 
Skyrme and Gogny effective nuclear interactions and focus on the possible existence of 
permanent imanation in the neutron plasma at finite temperature. The results show that, in 
the low density region, it is possible to have an increasingly stronger permanent polarization 
in the plasma at moderate temperatures as magnetic field strengths, higher than B w 10 16 G, 
exist. Below this strength the effects are almost negligible. At intermediate densities there 
is a plateau where both forces predict that the polarization remains approximately constant 
with maximum values around 12%. For higher densities, up to the limit where the nucleonic 
scenario holds, there is a dramatic difference in the behaviour predicted by these two forces. 
While Skyrme forces allow a ferromagnetic transition Gogny forces prevent it. The structure 
of the paper is as follows. In section [Til we describe the different nuclear potentials used 
in the description of the nucleon-nucleon (NN) interaction, namely the Skyrme and Gogny 
forces. Later in this section we discuss the thermodynamical observables involved in the 
description of a nuclear plasma in the presence of a magnetic field. Results will be presented 
in section Ifflt and summary and conclusions will be given in section IIVI 

II. THERMODYNAMICAL POTENTIALS AND NUCLEAR INTERACTION 
MODELS 

We will start by describing the system we want to study. We consider a pure homogeneous 
neutron plasma which is composed of particles with spin projection on the z-axis, a. It can 
be either a = +1 for spins aligned parallel to a uniform magnetic field that will be taken 
along the z-direction, B = Bk, or a = —1 for antiparallel spins. For a more complete 
description additional components in the hadronic and leptonic sector should be considered 
when beta equilibrium is imposed, however, for simplicity, we have restricted ourselves to 
the pure one component neutron plasma in this work. 

The neutron plasma density, p, is composed of up(+) and down(— ) spin aligned particles 

P = P+ + P- (1) 
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where 



Pa 



d 3 k 



Uk) 



(2) 



and the f a (k) is the Fermi distribution function for the fermion component with spin pro- 
jection a at temperature T = 1/(3 {k B = 1) 



fa(k) 



(3) 



X _|_ e /3{e a (k)-fi) 

The single particle energy e a (k) will be especified for each interaction model that will be 
used in this work, see below. The magnetization in the neutron plasma is defined as 



where \i 



M = fi n Ap (4) 
-1.9130427(5)/Zjv is the neutron magnetic moment in units of the nuclear 



magneton 21[ and A is the spin excess or polarization of the system 



(5) 



A. Skyrme force 

We have considered, in first place, the phenomenological Skyrme interaction that appears 
in the literature under a rather general form 22] 
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to (1 + x P a ) S(r) + ~ti (1 + xiP CT ) [k /2 5(r) + 5(r)k 2 

+t 2 (1 + x 2 P CT ) k' ■ 5(r)k + ~t 3 (1 + x 3 P CT ) P°(R)5(r) 
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(6) 



where r = ri — r 2 and R = (ri + r 2 )/2, k = (Vi — V 2 )/2z the relative momentum acting 
on the right and k' its conjugate acting on the left. P a is the spin exchange operator. Note 
that we have omitted terms not relevant for homogeneous systems. As widely known, this 
potential allows for a good reproduction of finite nuclei and their excited states and also 



bulk matter relevant to neutron stars 
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In order to study the neutron plasma we must consider the partition function, Z, of the 
available configuration space, Q. The one that must be used for our system with hamiltonian 
7i is that of the grand canonical ensamble under the influence of an external magnetic field. 



-PF M (n) 



Tr{e 



-pH(n) 



} 



(7) 
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The thermodynamical potential relevant for the description of a plasma at temperature 
T and external magnetic field H is the Helmholtz free energy, Fm, defined as 271 ] 



E-TS- ME (8) 



Note that previously we have used B to designate the magnetic field strength, however 
one should have in mind that apart from the external currents providing the external field 
H, the magnetization may contribute to the internal magnetic field [28] 

B = H + AnM (9) 

A self-consistent calculation would require considering the contribution of the magnetization 
to the total magnetic field B, however we expect M to be much less than the external field, 



so we will consider that the ratio | H/B | will always be close to unity [12j. In what follows 
we will keep the notation using B to designate the total magnetic field strength. In addition 
we will use capital letters for any extensive magnitude, A, and, for its associated intensive 
magnitude (per particle), the same letter in lower A/N. 

Using the non-relativistic Brueckner-Hartree-Fock approximation, the single particle en- 
ergy for a neutron with momentum k and spin projection a in presence of a magnetic field 
B can be written as 

e a {h) = %^ + U a (k)+n n aB (10) 
2m 



where U a (k) is the single particle potential For the Skyrme interaction it has a 

quadratic k-dependence that can be included into an effective mass plus a k independent 
term, U a , as 

h 2 k 2 — 

e a = —- + U a + fi n aB (11) 
2m* 

where the effective mass is given explicitly by 

m* 1 



- Xx) + t 2 {l + X 2 ))p-a + 2t 2 (l + X 2 ) Pa 



(12) 



Wr 

that in this approach is independent of temperature. It is useful at this point, in order to 
calculate the thermodynamical quantities of interest in the magnetized plasma, to define an 
effective neutron chemical potential into the following form 

Vefa = H-U a -n n aB (13) 
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Then, the Fermi distribution reads like that of a free Fermi gas with effective mass m* 
and chemical potential u e f a . 



m = — * - - (14) 

1 + e^a^r-^/') 

However, the single particle potential will be in general a function of k, U a (k), as i.e. 
when using the Gogny force. See later in this section. 

In the plasma the total energy per particle, for a given interaction model, can be obtained 
from averaging the kinetic, potential and field contributions 

1 r rl 3 h h 2 h 2 1 r rl 3 krl 3 h' R 2 

e=-E / iV-^Uk) + - E / -i^(V NN ) A Uk)Uk>) + f- (15) 
P V 7 ( 27r ) 2m P^ ( 27r ) 8?r P 

For the Skyrme interaction it can be casted into the following form: 

Skyrme 1 _ , ^ ^ _ , 

^+ + — r^- + 



p 2m!j„ p 2m*_ 



1 

P 

where r a can be expressed as 



t {l-X ) + ];t 3 {l-X 3 )p a 
O 



2 



^ + ^ (16) 



V/,(/.'! (17) 
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We have considered in this work two Skyrme parametrizations given by the Lyon group, 



namely we have chosen the SLy4 and SLy7 parametrizations [23|, [2J] as they are the most 
commonly used in the literature. These provide good values for binding of nuclei and neutron 
matter equation of state (EOS) giving values of maximum neutron star masses around 
1.5M & . In Table [J we summarize the values of some observables: saturation density, po, 
binding energy for symmetric nuclear matter, a v , symmetry energy, a s , and incompressibility 
modulus, K^, for the effective interaction models used in this work. 



B. Gogny force 



The Gogny force has been extensively used in the literature not only for describing finite 



nuclei but for the EOS of pure neutron matter 



31 



32 



331 ] . The interaction potential is 



given by the contribution of finite range terms and zero-range terms: 



V, 



Gogny . 



NN 



ri,r 2 ) = E{[W* + B t P° - H t P T - M t P° P T ] e -^-^ 2/ ^ + 



t 3l (l + x 3 P CT )p Q ^(r)} 



(18) 



TABLI 
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I: Values of some observables with the Skyrme and Gogny forces considered in this work 
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Model 


Po{fm~ 


" 3 ) (MeV) 


a v (MeV) 


a s (MeV) 


SLy4 


0.160 


230.9 


-15.97 


32.00 


SLy7 


0.158 


229.7 


-15.89 


31.99 


D1S 


0.1625 


203 


-16.01 


31.13 


DIP 


0.1737 


266 


-16.19 


34.09 



where P a and P T are the spin and isospin exchange operators. The values of some observ- 
ables for the two parametrizations considered in this work for the Gogny force appear in 

Table I. The D1S parametrization was orig inally introduced to describe the pairing proper- 

I 

ties and surface effects in finite nuclei [34J while the DIP aims at reproducing the EOS of 



neutron matter obtained using realistic interactions [35]. 

In this model the expression for the single particle potential felt by a neutron of momen- 
tum k and spin projection a is given by 

2 2 

U°° 9ny {k) = ^{t 3i (l-x )p ai p^ + a i t 3i p^ 1 } + 7r 3 / 2 J2 f if{[W i -H i ]p + 

i=l i=l 

[B t - M,K - [Wi +Bi-Hi- Mi]I[(k, T) - [Bt - M i ]I~' T (k, T)} (19) 

where the value of the If integral is 

i nk,T)=J^Uk')e-^ 2 (20) 

The integral I~ a in Eq. (fl~9l must be calculated using the oppositely polarized particle 
distribution f_ a . For the Gogny force we must note that the single particle energy entering 
the distribution function f a is given by 

e a (k) = ^ + XJ G a °* n \k) + ix n oB (21) 



Recent works [36( have shown that nucleon effective masses in plasmas described with 
Gogny forces have a weak temperature dependence thus, as an approximation, the most 
physically relevant interactions may be considered to take place around the Fermi surfaces 
with momentum kp a . Then, we will consider here a quadratic approximation in k-space to 
the full single particle potential since it proves to be a good approximation at least for k 



up to values around k m kp a as shown in ref. 29|]. In this way the quadratic momentum 



approximation for the single particle potential can be written as 

/ 1 rlTJ Go 9 n y((k^\ 

U^ ogny (k) « U° 09ny {k = 0) + -- AV , k2 ( 22 ) 

\2fc dk J k=kFa 

and an effective mass can be defined as 

m*_f m dU^v^ Y 1 

m-[ 1 + h 2 k> dk' ) k , =k {2S) 

To obtain the thermodynamical variables of interest once temperature, density, magnetic 
field strength and imanation are fixed, we must solve selfconsistently the set of equations 
given by 

m\ _(.,m dU^ n y(p,k',T,B) Y 1 

m-[ 1 + H*k dk' ) kl __ kFa (24) 

P = P+ + P- (25) 
Ap = P+ - p_ (26) 

Above, Eq. (12^1) represents the mean field felt by the up (down) particles in the plasma 
while Eq. (I25p must be solved in order to fix neutron number conservation and Eq. (|2"6"|) 
relates to the imanation constraints in the plasma. When analyzing the case of the Skyrme 
force the expression of the effective mass is analytical and Eq. ff24l) must be substituted by 
Eq. (HU). 

Once this set of equations is solved, we must evaluate other thermodynamical quantities 
on the hyper-surface A(p, T, B). The energy per particle is given by 

e Go 9 ny = +r _ ) + 

£ - x )p a '- 1 p +P - + ^ 2 rfm P ^±^[B t - Mi] - 

i=l P A Ap 

[W t + Bi — Hi — M 4 ]^[Gf' +) (p,T) + G { r>-\p,T)] - 
[Bt ~ Mi] 1 [G\ + '- ] ( P , T) + GY +) (p, T)] } + ^ (27) 
where the value of the Gi integral is 

G^'\p,T) = J ^- 3 U(k,T)lf(k,T) (28) 
The entropy per particle can be written then as 



s = — E t^/^)M/.W)-(i-/.(«i-/^)) (29) 

To extract the energetically favoured spin polarization in this ensamble the free energy 
per particle, /m, has to reach a minimum, A min , on the hypersurface at a given density, 
temperature and magnetic field strength. 

( df M \ 



OA 



I = (30) 

p,T,B 



III. RESULTS 



In this section we present the results. We have considered magnetic fields up to a max- 
imum value according to the scalar virial theorem. For a neutron star of radius R the 
maximum allowed stored magnetic energy is given by the relation ^^- B ^ x ~ ^jf- then, for 
a typical neutron star with M = 1.5M and R = 10~ 5 R Q this yields an interior maximum 
value of the magnetic field B max « 10 18 G. Recent observations report surface values up 



strengths have been speculated [371 ] . 



to B m 10 15 G as deduced from hydrogen spectral lines [3| but even higher magnetic field 



Skyrme models show a pathology, a ferromagnetic transition in the range [1 — 4]po, that is 
due to the appearance of a spontaneous polarization in the system even if no magnetic fields 
are present. To illustrathe this we show in Fig. [1] the development of a non-zero value of 
polarization A where the free energy per particle achieves a minimum until a ferromagnetic 
transition occurs around ~ 4p for the SLy4 model at a temperature T = 2 MeV and 
vanishing magnetic field strength. 

In Fig. [2] we show the free energy per particle for a plasma at p = 0.5po fc> r a typical 
temperature in a proto-neutron star, T = 10 MeV, and three different values of the magnetic 
field strength as calculated with the SLy4 parametrization. We can see that for fields below 
B m 10 17 G the effects introduced by the presence of a magnetic field are almost negligible. 
Then, when including the magnetic field the polarization symmetry is broken and no longer 
appear two minima. 

In the low density region the free energy is an increasing function of the density but the 
minimun of the polarization shifts to higher values as can be seen in Fig. [3] for T = 10 
MeV and B = 10 18 G using the SLy4 model. At densities below 0.5po a wide range of 
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FIG. 1: Free energy per particle in absence of magnetic field for different values of density and 
T = 2 MeV with the Skyrme SLy4 model as a function of the polarization. 
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FIG. 2: Free energy per particle at a density p = 0.5/9o for different values of magnetic field at 
T = 10 MeV with the Skyrme SLy4 model as a function of the polarization. 



models show the development of a frustrated state of matter that is known with the name 



of pasta ll7|,ll8|,ll9|,Q. Different simulations performed without the presence of a magnetic 



field show a variety of spatial shapes that form in neutron rich matter. It is important 
to explore the thermodynamical conditions for magnetized pasta where the presence of 
magnetic domains could largely impact neutrino opacities in the proto-neutron star early 
cooling stage 



In Fig. H]we show the polarization, in percentage (%), that minimizes the free energy at 
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FIG. 3: Free energy per particle for different values of the density at T = 10 MeV with the Skyrme 
SLy4 model and B = 10 18 G as a function of the polarization. 



T = 5 MeV resulting from a calculation with the SLy7 model as a function of density for p < 
4p - For higher densities a nucleon-quark transition is allowed by phenomenological models, 



see for example 



39(, so we will restrict our non-relativistic analysis to this range. We can 



see that there are three differentiated polarization regions. At densities increasingly smaller 
below ~ 0.5po ; the polarization rises with increasing magnetic field strength. It is worth 
mentioning that we have assumed that even at these low densities the gas picture still holds. 
However, the formation of fully spin saturated microscopic nuclear regions or "domains" 
S non>d be e.pncit.ye^ed^efutu.eH.Fo, densities in the range 0.5p ~ P ~ 3p a 
stable permanent imanation of the neutron plasma is allowed, if a magnetic field is present, 
giving a maximun polarization around 12% for the maximun B considered in this work. For 
densities larger than 3po ; again, there is a very rapid increase of the spin polarization that 
achieves full saturation and therefore a ferromagnetic transition occurs at around 4po for 
this particular parametrization. 

The robustness of the findings beforementioned can be seen on Fig. [5] where a comparison 
of Skyrme (SLy4 and Sly 7) and Gogny models (D1S and DIP) is shown for the case with 
p = po at T = 10 MeV and B = 10 18 G. We can see that all of them predict the same 
value of the minimun around A m j„ m 12%. The detail of the robustness at low densities 
with respect to the Skyrme (SLy7) and Gogny (DIP) parametrizations is shown in Fig. [6] 
for a density p = 0.5po at a T = 5 MeV where the permanent polarization is plot versus the 
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FIG. 4: Permanent polarization for different values of the magnetic field strength at T = 5 MeV 
with the Skyrme SLy7 model as a function of density. 

logarithm of B. Below B = 10 16 G the effect is negligible but it largely increases up to a 
value of ~ 12% for the B max . We can see that the increase of magnetization for sufficiently 
high magnetic field grows more rapidly than in the very low density and moderate magnetic 
field strengths regime where Curie's linear law M ~ B/T is expected to hold. 
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FIG. 5: Free energy per particle at a density p = po, T = 10 MeV and B = 10 18 G with different 
Skyrme and Gogny models as a function of the polarization. 

In Fig. [7] we show the plasma polarization at T = 5 MeV and B = 5 10 17 G with the 
Skyrme SLy7 (solid line) and Gogny DIP (dashed line) models as a function of the density. 
It is clearly observed the general feature of enhanced polarization at low densities and a 
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FIG. 6: Polarization at a density p = 0.5po an d T = 5 MeV with the Skyrme SLy7 and Gogny 
DIP model as a function of the logarithm of the magnetic field strength. 
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FIG. 7: Permanent polarization at T = 5 MeV and B = 5 10 17 G with the Skyrme SLy7 and 
Gogny DIP model function of the density. 



plateau at intermediate densities. However at high densities there is a dramatic difference 
since in the Skyrme model there is a ferromagnetic transition at ~ 4p while it is forbidden 
with the DIP model in agreement with realistic simulations [16] . In the case of the Gogny 
force there is a plateau with permanent residual polarization of about 2%. These results 
agree with those previously obtained with calculations performed in a relativistic scenario in 
the context of a mean field approximation where they estimate a maximum magnetization 
smaller than 4% 12]. 
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The effects of temperature on other thermodynamical quantities of interest can be seen in 
Fig. [HI Fig.[9]and Fig.[10]where the energy, entropy and free energy per particle, respectively, 
are shown for a density p = po ai B = 10 17 G at several temperatures for the D1S Gogny 
force. The typical range of temperatures relevant to proto-neutron stars goes up to values 
around T = 40 MeV. We can see that there is a moderate dependence of the energy per 
particle, e, with temperature. In the case of the entropy we can see the familiar decreasing 
behaviour (4l| with increasing order of the system, having a maximum around vanishing 
polarization. At high temperatures where correlations are greatly diminished one should 
recover the Sackur-Tretode formula limit s ~ 2.5 41]. 
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FIG. 8: Energy per particle at a density p = po for different values of temperature and B = 10 17 
G with the Gogny D1S model as a function of the polarization. 



In Fig. [IT] we show the entropy per particle as a function of the polarization for two 
different densities. The solid line corresponds to a density p = 0.5po and the dashed line to 
p = po a t T = 10 MeV with the Gogny DIP model. For this plot we have taken the upper 
limiting case B = 10 18 G. We can see that as density decreases in-medium effects become 
less important and allow the system more available configurations thus entropy increases. 

The effective mass for up and down polarized neutrons is shown in Fig. [12] for densities 
p = 0.5p (solid line) and p = p (dashed line) at T = 10 MeV with the Gogny DIP model. 
As polarization increases, the ascendent (descendent) curves for each density correspond to 
up (down) polarized neutrons. As shown, in-medium effects decrease the effective neutron 
mass as density increases. In Fig.[13]we show the effective neutron masses at a density p = po 
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FIG. 9: Entropy per particle, same thermodynamical conditions as in Fig. [HJ 



150 



100 - 



> 

CD 



50 



-50 







D1S 



-0.5 



T=0 
T=10 MeV 
T=20 MeV 
T=40 MeV 





A 



0.5 



FIG. 10: Free energy per particle, same conditions as in Fig. [HJ 



with the D1S model (solid line) and with the DIP model (dashed line) at T = 10 MeV as 
a function of the polarization. We can see that the effects due to the parametrization or in- 
medium correlations are much more important than those related to the polarization since 
for the maximum B considered in this work induces a A m j n « 12% that relates to a mass 
shift of an amount of 5m* « 5 MeV with respec to the unpolarized case, and 5m* ~ 10 MeV 
for the p = 0.5po case. 

The composition of the polarized plasma can be obtained by solving the set of eqs. ()24l) . 
(J25]) and (|26|) . In Fig. [HJwe show the polarized effective neutron chemical potential at densi- 
ties p = 0.5po (solid line) and p = po (dashed line) at T = 10 MeV and B = 10 17 G with the 
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FIG. 11: Entropy per particle for p = 0.5po (solid line) and p = po (dashed line) at B = 10 18 G 
and T = 10 MeV with the Gogny DIP model. 

1000 

950 

_ 900 
> 

CD 

850 

* 

E 800 
750 
700 

-1 -0.5 0.5 1 

A 

FIG. 12: Neutron effective mass for p = 0.5po (solid line) and p = po (dashed line) at T = 10 
MeV with the Gogny DIP model. The ascendent (descedent) curves are for up (down) polarized 
particles. 

Gogny DIP model as a function of the polarization. For each density as polarization grows 
the ascendent (descendent) lines show the up (down) polarized particle effective chemical 
potential. As the system is populated with more particles aligned (antialigned) with the 
magnetic field the energetic cost of aligning a new particle is bigger (smaller). As density 
increases this cost is also bigger due to nuclear correlations as extracted from the figure. 
In order to compare Skyrme and Gogny polarized population yields, in Fig. [TS] we show 
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FIG. 13: Effective neutron mass at a density p = po with the D1S model (solid line) and with the 
DIP model (dashed line) at T = 10 MeV as a function of the polarization. The ascendent (descent) 
curves are for up (down) polarized particles. 
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FIG. 14: Polarized effective neutron chemical potential at densities p = 0.5po (solid line) and 
p = po (dashed line) at T = 10 MeV and B = 10 17 G with the Gogny DIP model as a function of 
the polarization. See text for details. 



the polarized effective chemical potential at a density p = 0.5po for the SLy7 (solid line) 
and DIP (dashed line) models with B = 5 10 17 G as a function of the temperature. Up 
(down) effective neutron chemical potentials are depicted for each model as the upper (lower) 
curves. As the system is heated the permanent polarization stays rather constant however 
the correlations are diminished and the overall energetic cost is smaller. 
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FIG. 15: Polarized effective chemical potential at a density p = 0.5po with the SLy7 model (solid 
line) and the DIP model (dashed line) and B = 5 10 17 G as a function of the temperature. 

In Fig. [16] we show the polarized effective chemical potential densities p = 0.5po (upper 
panel) and p = Po (lower panel) with the SLy7 model (solid line) and the DIP model (dashed 
line) at T = 20 MeV as a function of the logarithm of the magnetic field strength. For each 
calculation model upper (lower) curves correspond to up (down) polarized particles. We can 
see how the effect of magnetic field strengths smaller than B pa 10 16 play almost no role in 
the plasma polarized neutron population. As B increases there is a gradual enhancement 
(decrease) of the up (down) polarized neutron sector. However the parametrization used 
plays an important role and Skyrme models predict systematically bigger effective neutron 
chemical potentials than the Gogny forces. From the comparison of upper and lower panels 
we can see that in-medium effects largely increase the effective chemical potentials. 

IV. SUMMARY AND CONCLUSIONS 

We have studied the effect of a strong magnetic field on some thermodynamical observ- 
ables in a pure neutron gas within the framework of the non-relativistic Brueckner-Hartree- 
Fock approximation. We have considered magnetic field strengths up to the maximum value 
energetically allowed by the scalar virial theorem B max pa 10 18 G. We have found that for 
magnetic field strength values smaller than B pa 10 16 G and due to the characteristic tiny 
value of nucleon magnetic moment the effect on plasma polarization is almost negligible. We 
have performed a comparative study of the plasma observables using effective Skyrme and 
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FIG. 16: Polarized effective chemical potential at 
densities p = 0.5po (upper panel) and p = po 
(lower panel) with the SLy7 model (solid line) and 
the DIP model (dashed line) at T = 20 MeV as 
a function of the logarithm of the magnetic field 
strength. For each model upper (lower) curves 
correspond to up (down) polarized particles. 

Gogny nuclear interaction forces in the density, temperature and magnetic field strength pa- 
rameter space. With these effective interaction models we find that in the low density limit 
increasingly large spin polarizations are allowed as density decreases from an upper bound 
of w 0.5p robustly for both models. Then, for intermediate densities, there is a rather 
constant polarization plateau where a permanent imanation saturates to a value of up to 
12% for the maximum B considered in this work and moderate temperatures. In the high 
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density region, up to where the meaningful hadronic picture holds, Skyrme and Gogny forces 
behave very differently. While the Skyrme interaction allows a ferromagnetic transition in 
the high density limit, the Gogny force prevents it. Temperature effects on magnetization 
remain moderate in the explored range of temperatures, up to about T = 40 MeV rele- 
vant to proto-neutron stars. The influence of a strong magnetic field should be considered 
increasingly relevant as density decreases, that is, in the description of the outer shells of 
proto-neutron stars below ~ 0.5po where a more detailed analysis must be performed in 
order to see ther role played by frustration. For bulk matter in the intermediate and high 
density range a relatively small polarization is allowed that only affects moderately the ther- 
modynamics. Beta equilibrium should also be appropiately treated since the appearance of 
electrically charged particles may affect the onset densities of enhanced polarization in the 
plasma. 
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